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Abstract 

Spatially extended dynamical systems, namely coupled map lattices, driven 
by additive spatio-temporal noise are shown to exhibit stochastic synchro- 
nization. In analogy with low-dymensional systems, synchronization can be 
achieved only if the maximum Lyapunov exponent becomes negative for suffi- 
ciently large noise amplitude. Moreover, noise can suppress also the non-linear 
mechanism of information propagation, that may be present in the spatially 
extended system. A first example of phase transition is observed when both 
the linear and the non-linear mechanisms of information production disappear 
at the same critical value of the noise amplitude. The corresponding critical 
properties can be hardly identified numerically, but some general argument 
suggests that they could be ascribed to the Kardar-Parisi-Zhang universality 
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class. Conversely, when the non-linear mechanism prevails on the linear one, 
another type of phase transition to stochastic synchronization occurs. This 
one is shown to belong to the universality class of directed percolation. 
PACS numbers: 05.45.-a,05.45.Jn,05.45.Ra,05.45.Xt,05.40.Ca 
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Introduction 



The synchronization of dynamically coupled chaotic systems is a new interesting phe- 
nomenon that has attracted many studies and efforts in the last years. The first models 
and examples discussed in the literature concerned the case where the dynamical coupling 
is performed by a deterministic chaotic signal (either external or self-generated) acting on 
different trajectories of the same low- dimensional chaotic system (e.g., see Under quite 
general conditions, e.g. large enough coupling strength and sufficiently long evolution time, 
the two trajectories approach each other and, despite their chaotic nature, no more diverge. 
Beyond the intrinsic conceptual interest of " chaotic synchronization" , the main motivations 
of these studies came from the possibility of identifying efficient mechanisms for controlling 
or suppressing chaos in deterministic signals, with promising perspectives of applications to 
the automatic control of devices and to cryptography. 

On the other hand, applying a control mechanism to a deterministic system by a chaotic 
signal looks artificial with respect to the more realistic scenario, where the coupling could be 
produced by noise. It is worth stressing that deterministic chaos cannot be assimilated to a 
random process, despite in some cases they may yield equivalent results []. Accordingly, it is 
quite natural to raise the question if synchronization can be observed also when the chaotic 
coupling is substituted with a stochastic signal. The problem of "stochastic synchronization" 
has been investigated by several authors considering two different initial conditions of 
the same low-dimensional chaotic system (e.g. the logistic map or the Lorenz system) and 
making them evolve by adding to both of them the same realization of noise. For sufficiently 
large values of the noise amplitude, the stationary probability measure associated with 
the chaotic dynamics may be modified by the addition of noise favouring the stabilization 
onto the same stochastic orbit. As argued by Pikovsky ||, this effect can be explained 
and quantitatively predicted by observing that, at some value of the noise amplitude, the 
Lyapunov exponent associated with the dynamics changes from positive to negative. A 
necessary condition for the occurrence of synchronization is that the dynamical system has 
expanding and contracting regions, so that noise may amplify the role of the latter against 
the former. For instance, synchronization cannot emerge from uniformly hyperbolic maps, 
like the Bernoulli shift. In Section I we discuss "stochastic synchronization" in simple maps, 
that we introduce as examples of different classes of models. The new question that we want 
to address in this paper is the following: under which conditions stochastic synchronization 
may occur in spatially extended dynamical systems coupled by the same realization of spatio- 
temporal noise? The models, that we shall study numerically, are the diffusively coupled map 
lattice (CML) [^] versions of those introduced in Section I. A comparison with the analysis 
therein worked out on stochastic synchronization in low-dimensional systems shows that the 
spatial structure makes this phenomenon even richer. Specifically, let us list hereafter the 
main results contained in this paper: 

• Pikovsky's criterion can be generalized to spatially extended systems : Stochastic 
synchronization occurs only if the maximum Lyapunov exponent A M is negative. In 



For instance, the Monte Carlo procedure applied to an ergodic system yields the same equilibrium 
properties emerging from the Liouville measure associated with the dynamics. 
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fact, we find that for noise-coupled CML models A becomes negative for values of the 
noise amplitude a larger than a threshold value <ta- 

It is well known that in spatially extended systems information may propagate with 
a finite velocity V 0, due to spatio-temporal instabilities. For noise-coupled CML 
models, with smooth nonlinearities, the only relevant instability mechanism is the 
linear one, so that V vanishes with A. This identifies a synchronization transition 
(PT1) at a threshold value oy = a\, whose critical properties seem to belong to the 
Kardar-Parisi- Zhang (KPZ) universality class of interface roughening transition [TT)|. 



We discuss also other models where deterministic dynamics or spatiotemporal noise 
may induce strong non-linear effects. In these cases we find that V vanishes at a value 
of the noise amplitude ay > o~a- A different kind of synchronization transition (PT2) 
is observed at ay and its critical properties are found to belong to the universality 
class of directed percolation (DP) fllTt|l2 . 



In the latter case we find also that the time needed to achieve a synchronized state 
grows exponentially with the system size L for ay > a > a\ : this dynamical regime is 



analogous to the "stable-chaotic" phase observed in Ul3| . Conversely, this time grows 
logarithmically with L for a > ay (see also 

In Section II we introduce the models and the suitable indicators for describing the syn- 
chronization transitions induced by spatio-temporal noise in CMLs. The dynamical features 
of the different scenarios are analyzed in Section III. In Section IV, we study the critical 
properties associated with the transition to stochastic synchronization in spatially extended 
systems : the numerical analysis indicates that this phenomenon is a genuine nonequilibrium 
process, characterized by the scaling laws of KPZ and DP universality classes for PT1 and 
PT2, respectively. We expect that a suitable coarse-graining or renormalization approach 
could provide a rigourous ground for these results if one could map the noise-coupled CML 
dynamics onto the stochastic PDE's that identify the above mentioned universality classes. 
Unfortunately, we have not been able to work out any consistent derivation. Nonetheless, 
numerics allows one for drawing some conclusions about the problem at hand: they are 
summarized in Section V, together with some perspectives for future work on this subject. 



I. A SURVEY ON STOCHASTIC SYNCHRONIZATION IN 
LOW-DIMENSIONAL SYSTEMS 

Before entering the main topic of this paper, i.e. stochastic synchronization in spatially 
extended systems, it is worth presenting concepts, tools and properties characterizing the 
same kind of phenomenon in low-dimensional dynamical systems. In this section we consider 
some simple models, that exemplify the general features of this phenomenon. 

The basic model equation is the stochastic map 

x^ 1 = + a ■ rf (1) 

where the state variable x t is a real quantity depending on the discrete time index t, a is the 
amplitude of the time dependent random variable 77*, uniformly distributed in the interval [- 
1,1], and f(x) is a map from S C M into the interval I C M. Stochastic synchronization can 
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be investigated by considering two different initial conditions, x° and y° of dynamics ([!]) with 
the same realization of additive noise rf. More precisely, we assume that the corresponding 
trajectories, x l and y*, generated by (§) may synchronize if their distance 

d(t) = \x t -y t \ (2) 

becomes smaller than a given threshold A, usually assumed much smaller than unit. This 
definition allows one to identify two natural quantities associated with the synchronization 
of trajectories: the first passage time Ti(A), i.e. the first instant of time at which d(t) < A 
and the synchronization time t 2 (A), i.e. the interval of time during which d(t) remains 
smaller than A |f[5|| . We want to point out that, in general, both Ti(A) and r 2 (A) depend 
on the initial conditions and on the realization of noise. Accordingly, their averages over 
both ensembles have to be considered as the quantity of interest for our analysis. For the 
sake of simplicity we shall indicate these averaged times with the same notations Ti(A) and 
r 2 (A). 

In all the numerical simulations, for sufficiently small value of A (typically A < 10 7 ), 
the results do not depend on the choice of A. Moreover, one can assume that the trajectories 
have eventually reached the synchronized state if, after Ti(A), d(t) remains smaller than A 
for arbitrarily large integration times. This amounts to say that t 2 (A) goes to infinity with 
the integration time. This heuristic definition of the synchronized state can be replaced 
by a more quantitative criterion: according to Pikovsky 0, for a given value of a the 
synchronized state occurs if the Lyapunov exponent A of dynamics (|]) is negative. This 
indicator is defined as follows: 

A = lim -InT] \ JLr\ (3) 
where the dynamical variable obeys the "linearized" dynamics 

and the derivative is computed along the trajectory given by ([!]) . We want to remark that 
A is well defined for deterministic dynamics (e.g. the a = case in (|Tj) ). In the framework 
of the linear stability analysis a positive (negative) A measures the average exponential 
expansion (contraction) rate of nearby trajectories. For a/Owe are faced with stochastic 
trajectories and it is not a priori obvious if A is still a meaningful quantity ||16|| . On the 
other hand, eq.(§) does not depend explicitely on 77* and it is formally equivalent for both 
the noisy and the noise-free cases. However, the presence of noise modifies the evolution of 
the system w.r.t the noise-free case and, accordingly, also the tangent space dynamics. We 
have verified numerically that A is a self-averaging asymptotic quantity also for dynamics 
(H|) with a 0. It can be interpreted as the average exponential expansion (contraction) 
rate of infinitesimal perturbations of stochastic trajectories generated by the evolution rule 
(H|) . In particular, its value is found to depend on a, but not on the realization of noise. 

Let us remark that it is quite simple to argue why A < implies r 2 — > 00: After some 
finite time Ti(A), d(t) has decreased below a small threshold A and the trajectories can be 
viewed as a perturbation of each other. Linear stability implies that their distance will keep 
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on decreasing exponentially with an average rate A, so that, within numerical precision, 
they will converge rapidly onto the same trajectory. 

As a first example, we consider the continuous map shown in Fig. [I]: 

{— ctanh(6 (1 + x)) if x < — 1; 
ax (1 — \x\) if \x\ < 1; (5) 

ctanh(6 (1 — x)) if x > 1. 

where 5 =R and X = [—1, 1]. We choose the parameter values a = 4, and c = 0.5, so that 
(H) can be viewed as a sort of anti-symmetrized version of the logistic map at the Ulam 
point, taking values over the whole real axis. It can be easily shown that for o = and 
independently of b, A = In 2, i.e. map (|5]) is chaotic. Notice that the noise term of amplitude 
a extends X to the interval [— 1 — a, 1 + a] . We have verified numerically that , for any value of 
b and for a larger than a threshold value a\, A becomes negative and after some finite time 
T\ a synchronized state is eventually achieved. In particular, we find that cta is strongly 
dependent on b: for instance, we have obtained oa = 1.2 for b — 2 and a a = 0.019 for 
b = 1000. As recently found also by Lai and Zhou for a similar mapping, this result 
indicates that symmetric , i.e. zero-average, noise can yield stochastic synchronization. 

It is worth mentioning that some time ago Herzel and Freund || conjectured that stochas- 
tic synchronization can be achieved only if noise has a non-zero average. They were led to 
such a conclusion by studying stochastic synchronization for the case of a map / of the unit 
interval into itself, i.e. S = X = [0, 1] . In such a case the application of the stochastic 
evolution rule (|]) demands the adoption of some further recipe for maintaining the state 
variable x l inside the unit interval, when o ^ 0. For instance, one can choose the following 
reinjection rule x t+1 — > x t+1 + 1 ( x t+1 — > x t+1 — 1) if x t+1 < (x t+1 > 1). As discussed in 
0, any recipe of this kind yields an effective state-dependent noise, that does not preserve 
the original symmetry of the stochastic process rf, thus acquiring a non-zero average value. 
The above described example and the results obtained in || disprove their conjecture. 

Nonetheless, an interesting observation is contained in [fj: the stochastic evolution rule 
(0) for maps definite on a finite interval induces strong non-linear effects due to the discon- 
tinuities introduced into the dynamics by the state-dependent noise. This strong non-linear 
character of the dynamics is irrelevant for stochastic synchronization in low-dimensional 
systems; conversely, it reveals a crucial property for discriminating between different critical 
behaviours in high- dimensional systems, as we shall discuss in Section IV. 

It is also worth considering that even the presence of non-zero average noise does not 
necessarily guarantee stochastic synchronization. For instance, a counterexample is provided 
by considering dynamics ([I]) for the logistic map at the Ulam point: 

f{x) = 4x(l - x) (6) 

whose noise-free dynamics is mixing. The state-dependent noise modifies the probability 
measure of x l in such a way to increase the weight of the contracting regions of the map. 
On the other hand, A remains positive as well as T2 remains finite for any value of a 6 
[0, 1], despite they can be made so small and so large, respectively, to produce an apparent 
synchronization effect. As discussed in misleading results can be obtained in this case 
due to the finite computational precision of numerical simulations. 
As a final example, we consider the map 
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r(\_jbx if < x < 1/6 ,_, 

J{X) ~ \a + c{x- l/b) ifl/6<z<l ^ J 

Its dynamics converges to a stable periodic attractor: for b = 2.7, a = 0.07 and c = 0.1 this 
is a period-3 orbit with negative Lyapunov exponent, A = ln(cb 2 ) ~ —0.316 . Numerical 
analysis shows that the Lyapunov exponent of the stochastic evolution rule (P remains 
negative for any value of a and, according to Pikovsky's criterion, a synchronized state is 
always achieved, as we have numerically checked. 

In summary, we have presented in this section three different examples of low-dimensional 
dynamical systems epitomizing possible scenarios concerning the phenomenon of stochastic 
synchronization: In the last example it occurs for any value of the noise amplitude a, in the 
second one for no value of a and only above some threshold value, a\, in the first example. 



II. STOCHASTIC SYNCHRONIZATION IN CML MODELS 

The generalization of the stochastic map (JT|) to a CML model J7| with additive spatio- 
temporal noise can be defined by the following two-step evolution rule: 

x\ = {l-e)x\ + ~{x\_ x + x\ +1 ) 

x^ = m) + a- V t (8) 

The real state variable x\ depends now also on the discrete space index i = 1, 2, • • • , L: 
assuming unit lattice spacing, L corresponds to the lattice size. The strength of the spatial 
coupling between nearest neighbour maps in the lattice is fixed by the parameter e, that can 
take values in the interval [0,1]. Note that this kind of coupling amounts to a discrete version 
of a diffusive term. The application of the map / mimicks the reaction term of reaction- 
diffusion PDE's, so that CMLs are commonly assumed to represent a sort of discretized 
version of continuous PDE's. At variance with standard CML models, dynamics (§) contains 
also a stochastic term given by a set of identical, independent, equally distributed (i.i.e.d.) 
random variables {//«}, whose amplitude is determined by the parameter a. In what follows 
these i.i.e.d. random variables are assumed to be uniformly distributed in the interval [-1,1]. 

In full analogy with the low- dimensional case discussed in the previous section, stochastic 
synchronization can be investigated by considering two different initial conditions {x®} and 
{y®} for dynamics (§), coupled by the same realization of additive spatio-temporal noise. 
More precisely, we assume that the corresponding trajectories, {x\} and {?/*}, generated by 
(§) may synchronize if their distance 

m = ^Stf-ril (9) 

becomes smaller than a given threshold A, usually assumed much smaller than unit. Upon 
this definition one can straightforwardly extend to the CML case indicators like the first 
passage time Ti(A) and the synchronization time T2(A). Also in this case we assume that 
the trajectories synchronize if after the time Ti(A), t 2 (A) is found to diverge with the 
integration time. In all the numerical simulations we have used the same value of A as in 
the low- dimensional case. In fact, we have verified that also for the high-dimensional case 
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the results of numerical simulations are not affected by the choice of A, provided it is small 
enough, tipically A < 1(T 7 . A gain, both Ti(A) and T2(A) are quantities averaged over initial 
conditions and over realizations of noise. 

In analogy with low- dimensional systems, we expect that the suitable dynamical indi- 
cator for identifying stochastic synchronization in dynamics (^j) is the maximum Lyapunov 
exponent, A. In particular the phenomenon is expected to occur only for values of the noise 
amplitude a such that A is negative. For what concerns the interpretation of this indicator 
for the stochastic CML dynamics the same kind of remarks and conclusions discussed 
in Section I for the low dimensional case still hold. All the numerical estimates of A have 
been performed by applying the standard algorithm outlined in ||. 

Another relevant indicator, strictly related to the spatial structure of CML dynamics, is 
the average propagation velocity of finite amplitude perturbations || 

V = lim lim (10) 



where 



N{t) wit h in otherwise 



i=l 



is the number of non- synchronized or "infected" sites at time t. Here {xj} and {yj} represent 
the trajectories generated by dynamics @ starting from two initial conditions that differ of 
finite amounts 5i ~ 0(1) only inside a space region of size S: 

of i$ + 6 t if |L/2-i|<£/2; 
I x? otherwise 



The average < • > in ( fLOf) is performed over different initial conditions and noise realizations. 
The indicator V measures the rate of information propagation : we want to point out that 
V can take finite values even for non-chaotic evolution, i.e. for A < 0. For instance, 
this scenario has been observed for CMLs made of discontinuous maps [HJ. In this case 
infinitesimal perturbations are unable to be amplified, while finite amplitude perturbations, 
induced by the discontinuity, can propagate thanks to the spatial coupling. Therefore in 
spatially extended systems the information flow is absent only when V vanishes. 

As a first example let us consider the stochastic CML dynamics (g), equipped with 
map (|^) for b = 2.0 (see Fig.|l]). In this case, numerical simulations indicate that the 
stochastic synchronization of trajectories occurs for any value of the diffusive coupling e 
and for sufficiently large values of a, above which A becomes negative. For instance, with 
e — 1/3 we find that synchronization is obtained for a > a\ = 2.4768. This exemplifies the 
validity of Pikovsky's criterion also for spatially extended systems. 

Let us provide more details of the dynamics below and above the threshold value cta- 
For cr < cta we find that 7i(A) diverges exponentially with the system's size L, i.e. T\ ~ 
exp (L/£). The length scale factor £ is found to be independent of L and proportional to the 
inverse decay rate of the space correlation function of x\. Accordingly, L/£ is an estimate 
of the number of effective independent degrees of freedom. In this dynamical regime the 
probability P(£, A) that two trajectories get closer than a distance A is proportional to 
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the combined probability that each one of the L/£ degrees of freedom get closer than A, 
i.e. P(£, A) oc A L ^. One can reasonably assume rf oc P(£,A); this rough argument 
explains why T\ ~ exp(L/£). Note that even if d(t) eventually becomes smaller than A 
synchronization is rapidly lost because of the linear instability mechanism due to positive A 
and T2 is always finite. 

Conversely, a logarithmic dependence of Ti(A) on L is found for a > a a (see Fig. || (b)). 
Numerical simulations show that in this case the number of regions made of a few synchro- 
nized sites increases as time flows. Moreover, once a synchronized region is formed, it grows 
linearly in time, until all regions merge and synchronization sets in over the whole lattice. 
One can introduce an argument accounting for the logarithmic dependence T\ oc ln(L). An 
effective rate equation for the number of synchronized sites, n(t), can be constructed by 
assigning a probability p for the formation of new synchronized sites and a rate 7 for the 
linear increase of synchronized regions: 

Oti 

_ = 7+p (L_ n ) , (13) 

with < n(t) < L. This equation can be solved with the initial condition n(0) = , so that 
an estimate of T\ is obtained by imposing the condition n(ri) = L: 



n = — In 
p 



P - + l 

7 



(14) 



Note that a logarithmic dependence on L is consistent with the condition — >> 1, that can 
be satisfied for sufficiently large values of L. By estimating the parameters p and 7 directly 
from numerics we have verified that they fit reasonably well the simple phenomenological 
eq- (0). 

The indicator V does not provide any additional information about this kind of syn- 
chronization transition. Actually, V is found to be positive for a < <ta, while it vanishes at 
a = a a- This implies that in this model the linear mechanism of information production 
associated with a positive A rules also the propagation of finite amplitude perturbations |L7| . 



In summary, for a > o~a A is negative and V = , so that after the trajectories have reached 
7~i it seems that no information production mechanism can be responsible of the resurgence 
of d(t) above the threshold A. This is confirmed by numerical simulations, although in 
principle one could not exclude that the occurrence of a temporarily positive exponential 
expansion rate at some lattice site might produce a local amplification of d(t). In this respect, 
we want to remark that A is a global indicator, i.e. the exponential expansion rate between 
nearby orbits averaged in time and over the whole phase space. Accordingly a negative value 
of A is fully compatible with the above mentioned local, instantaneous event. This point, 
that reveals important for the understanding of the dynamical mechanisms underlying the 
synchronization transition, will be analyzed more carefully in the following section. 

Different scenarios are obtained by considering dynamics (|8|) with / given by the logistic 
map (@). For e = 1/3 and large enough values of L, one recovers features very similar to 
the case of a single logistic map, where the synchronization transition is absent. In fact, for 
any value of a, A and V are positive, while T\ and T2 remain finite. Moreover T\ is found to 
diverge exponentially with L with a parameter dependent rate. 



9 



A different situation occurs for e = 2/3, where A vanishes at a\ = 0.27, while V remains 
positive up to ay ~ 0.4 (see Fig. |3|). According to Pikovsky's criterion one expects that 
synchronization occurs above OA.This is indeed the case, although, for a A < a < a v , 7i(A) 
is still found to increase exponentially with L, while for a > ay , 7"i(A) grows logarithmically 
with L (see Fig. Despite the strong analogy between this transition at ay and the one 
occurring in the first example at a\, we want to remark that there is a crucial difference 
between them: below threshold r2 diverges in the former case, while it is finite in the latter. 

This indicates the existence of a new kind of synchronization transition at ay essentially 
ruled by V. Let us point out that, at variance with the first example, for a A < a < ay the 
non-linear mechanism of information propagation is enough for maintaining the exponential 
dependence of T\ on L. On the other hand, after t\ two trajectories get very close to each 
other and the negative A stabilizes them onto the same stochastic trajectory. 

Let us stress that the above described scenarios are not peculiar of the specific choice 
of the parameter values e = 1/3 and e = 2/3. We have checked that their main features 
are robust w.r.t. small but finite variations of s, although a systematic investigation of the 
parameter space would demand too huge computational times. 

Upon these examples one can conclude that in a CML of very large but finite size L 
stochastic synchronization of two trajectories within an accessible time span occurs when V 
vanishes. 

Moreover, a very similar situation is obtained when considering dynamics @ for the 
model of period-3 stable maps (0) introduced in |13[ . In this CML model A is found to be 



negative, independently of the value of the diffusive coupling parameter e and, accordinlgy, 
the dynamics eventually approaches a periodic attractor. On the other hand, the model 
exhibits a transition between a frozen disordered phase with V = to a chaotic phase with 
V > at e c ~ 0.6. The peculiar feature of this transition is that these two phases are 
separated by a small fuzzy region centered around e c , where both positive and null values 
of V can be observed up to available numerical resolution. 

The addition of noise to this CML dynamics according to eq. (§) has interesting conse- 
quences: A is kept negative independently of the noise amplitude a, while a small amplitude 
noise destabilizes the frozen disordered phase: V becomes positive even for e < e c , so that 
the fuzzy transition disappears. This notwithstanding, by increasing a up to a critical value 
ay(e) , V is found to drop again to zero not only below, but also above e c . For instance, 
one has ay ~ 0.16 for e = 0.58 and ay ~ 0.18 for e = 0.62. In both cases we recover the 
same kind of mechanisms characterizing the synchronization transition discussed for coupled 
logistic maps with e = 2/3. 

We want to point out that a finite value of V, when A is negative, is usually reported 



as a typical signature of a strong non-linear effect [l7j. For instance, it has been shown 
|i~3| that the discontinuity of map (0) yields such an effect already for the noise-free CML 
dynamics. Even if the discontinuity of map (0) is removed by interpolating the expanding 
and contracting regions with a sufficiently steep segment the effect is maintained |19| . 

The reinjection mechanism introduced by additive noise in maps of the interval [0,1] into 
itself produces similar discontinuities in the dynamics. As we have shown this non-linear 
effect is sufficient for giving raise to dynamical phases with A < and V > in the noisy 
dynamics (§). 

All these observations suggest to investigate if a similar scenario can be obtained for map 
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(§) in dynamics (§), by introducing an almost-discontinuity in the map, since in this case 
the reinjection mechanism is not present. This is easily obtained by taking a sufficiently 
large value of the parameter b, e.g. b ~ O(10 3 ). Numerical simuations show that, still for 
e = 1/3, it exists a range of a- values for which V > and A < 0. In Fig.[| we show the 
dependence of T\ on L in two regions of the parameter space where A is negative, while V 
is either positive or null. If V > an exponential increases of T\ with L is again observed, 
while a logarithmic dependence characterizes the dynamical phase with V = 0. 

Accordingly, this second kind of transition is not just an artifact due to discontinuities 
introduced by the rejection rule, but a mere consequence of sufficiently strong non-linear 
effects, that may be produced also in the absence of any rejection mechanism and also if the 
noise distribution maintains its symmetry. 



III. CRITICAL PROPERTIES OF THE SYNCHRONIZATION TRANSITION IN 

CMLS 

In the previous Section we have described two different kinds of phase transitions for 
stochastic synchronization where the noise amplitude a plays the role of the control pa- 
rameter. When both A and V vanish at a = cta the transition is from a non- synchronous 
dynamical phase to a synchronous one: we have denoted it with PT1. In the other case, that 
we indicate with PT2, A passes from positive to negative values while V remains positive up 
to a = ay, where we have observed a transition between different synchronous dynamical 
phases, characterized by an exponential and by a logarithmic dependence on the system size 
of the first passage time T\ . 

In both cases we have found that the correlation length £, defined in Sect. II, diverges 
at the transition point : this is an indication in favour of a continuous phase transition. 
Nonetheless, l/£ is found to vanish less than linearly when a — > cr^ in PT1, while it vanishes 
linearly in PT2 when a — > oy (see Figs. |6| (a),(b)). 

Moreover, close to the critical point the time averages of the mean distance between 
trajectories, z l = ~ Yh=i \ x \ ~ v\\i an d of the topological distance, p(t) = j Y,f=i h\ = 1 
if \zj\ > A, otherwise h\ = ) exhibits a continuous dependence on a (for the sake of space 
in Fig.^|we show these quantities only for the case of coupled logistic maps). 

We want to remark that our definition of the synchronized state implies that it exists a 
stable stochastic orbit that prevents the trajectories of the dynamical system to flow apart 
from each other below some very small, but finite threshold A so that r 2 = oo. In both 
cases numerical simulations show that this is what happens above the critical point, where 
A < and V = 0. If at some lattice site % a fluctuation makes the local Lyapunov multiplier 
positive (i.e., ln|/'(x*)| > 0, where /' indicates the first derivative of the map w.r.t. its 
argument), giving rise locally to the exponential divergence of nearby orbits, the process is 
rapidly reabsorbed due to the lack of any mechanism of information propagation. In this 
sense the synchronized state should be equivalent to the absorbing state, typical of Directed 



Percolation (DP) processes |TT|,[T2" 



Below the critical points the role of fluctuations determines the difference between PT1 
and PT2. The inspection of the space-time evolution of dynamics (H), using the symbolic 
representation of the state variables hj, is quite helpful for visualizing such a difference. For 
what concerns PT1, when a — > one observes that non-synchronized clusters propagate 
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as time flows (since V is positive) and, even if some of them may eventually die, in the 
meanwhile new ones have started to propagate, emerging also from previously synchronized 
regions. Since also A > any local fluctuation of d(t) produced by a positive multiplier has 
a finite probability to be amplified and eventually propagated through the lattice. 

On the contrary, sufficiently close to PT2, for a — ► oy, non-synchronized clusters never 
emerge from already synchronized regions and any connected non-synchronous cluster even- 
tually dies at Ti(L, a) in a lattice of finite size L (a situation very similar to what is observed 
in the active phase of DP as a finite size effect). Even if, in principle, the non-linear mecha- 
nism of information propagation is active, this suggests that, inside an already synchronized 
region, any local fluctuation of the Lyapunov multiplier towards positive values never persits 
long enough to activate the non-linear process of information propagation. 

Upon these numerical observations and exploiting the analogies with other critical phe- 
nomena, we are led to conjecture not only the existence of an absorbing, i.e. synchronized, 
state but also that PT2, at variance with PT1, should belong to the universality class of 
DP. This can be confirmed only by direct measurements of the critical exponents associated 
with the synchronization transitions. It is worth stressing that numerical estimates have 
been performed by approaching the critical points from below, i.e. for V —>■ + . 

Here we report the analysis of PT2 in the case of coupled logistic maps with e — 2/3. As 
usual a reliable measurement of any critical exponent demands a very accurate estimate of 
the critical point, i.e. ay, in this case. For this reason we have performed careful simulations 
for evaluating the dependence on the system size L of t±, that corresponds to the absorption 
time in the DP language. At the critical point a = ay, this time should diverge as 

n(L,a v )~L* (15) 



where z = u\\/u± is the dynamical exponent [11.12]. The quantity Tx(L,o~) is reported in 



Fig.|S] as a function of L in a log-log scale for different values of a. The best scaling behaviour 
is obtained for ay = 0.4018 , where one has z = 1.55 ± 0.05. This result agrees quite well 
with the most accurate numerical estimates of the DP value z = 1.5807 ||20|| . Relying upon 



this result, we have also measured the critical exponent associated with the temporal decay 
of the density p(t) of active sites, i.e. those sites where \x\ — yj\ > A . The DP transition 
exhibits at the critical point the following scaling law for the topological distance: p{t) ~ t~ s , 
with 5 = /3/v\\ = 0.1595 [pOj] . The density p = p(t) is shown in Fig. |9] for three different 
system sizes, namely L = 500, 1500 and 2000. As expected the scaling region increases 
with L and for L = 2000 an optimal fitting in the interval 3.2 < log(t) < 4.8 provides the 
estimate for the critical exponent S = 0.159 ±0.002. This confirms that PT2 belongs to the 
universality class of DP. 

For what concerns PT1 we have considered coupled maps of the type (^]) for b = 2.0 
and e = 1/3. The best scaling for the T\ as a function of L (accordingly to equation (|T5|)) 
is observed for a noise amplitude <j\ = 2.5015. For a = <j\ we have obtained the following 
estimates for the critical exponents : z ~ 1.01 — 1.04 and 5 ~ 0.35. Such values certainly do 
not correspond either to DP, or to any known universality class for percolation or growth 
processes. This seems analogous to what has been pointed out by Grassberger \T2\ for 



systems exhibiting incomplete deaths: when the asymptotic state is not a truly absorbing 
one the critical properties of DP cannot be recovered. 

Since in this situation an absorbing state seems not to exist, the dynamical exponent 
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z can be better estimated by measuring the number of "infected" sites of the chain N(t) 
defined in ([H]). At the critical point a = the following scaling law is expected to hold 

N(t) ~ t 1/z . (16) 



The data from numerical analysis are reported in Fig. [10]. For short times (t < 500) we 
obtain the inverse of the the dynamical exponent 1/z ~ 0.47 ± 0.05, a value consistent with 
the one expected for the Edwards- Wilkinson universality class zew — 2.0 | ]21|| . For longer 
times we observe a crossover to a lower z-value that is consistent with the one expected for 
the Id KPZ universality class (namely zkpz = 3/2). We think that these results could be 



interpreted in terms of the conjecture reported in [p3fl . In that paper it has been suggested 
that, for generic synchronization transition of coupled spatio-temporal chaotic systems with 
continuous state variables, the appropriate universality class should be the one of the KPZ 
model with a non-linear growth limiting term This idea originates from the obser- 

vation that, close to the synchronization transition, it is possible to describe the dynamics 
of small perturbations in terms of a reaction diffusion model with multiplicative noise |12| . 



Finally, this model can be mapped via a Hopf-Cole transformation into an equation that 
corresponds to a KPZ model with a non-linear term that prevents the surface from growing 



indefinitely . The critical scaling laws for this kind of models have been reported in |22 
the dynamical exponent z is found to coincide with the KPZ-one, while the other expo- 
nents are found to be different from the standard KPZ ones. Unfortunately a full numerical 



comparison of PT2 with the model discussed in |22[] is prevented by two major technical 
problems: the extreme difficulty to estimate with sufficient precision the critical value <ta 
and the strong finite size effects. 



IV. CONCLUDING REMARKS 

We have studied the problem of stochastic synchronization induced by additive spatio- 
temporal noise in CML models. In analogy with the low- dimensional case, synchronization 
of trajectories is observed if the maximum Lyapunov exponent becomes negative at a critical 
value of the noise amplitude. We have also identified two different critical behaviours asso- 
ciated with the synchronization transition. One of them belongs to the universality class of 
Directed Percolation. While the other one is not clearly identified, but indications suggest 
that it could belong to the universality class of the KPZ model with a non-linear growth 
limiting term. 

In particular, the DP-like phase transition describes the crossover between an "active" 
and an "absorbing" phase. The former is characterized by an exponential dependence on 
the system size L of the time needed for achieving the synchronized state, while the latter 
exhibits a logarithmic dependence. This scenario is reminescent of the phenomenology asso- 
ciated with stable-chaos in CMLs ||13|| , where the dynamics approaches a periodic attractor, 
rather than a "synchronized" stochastic trajectory. 

The two kinds of synchronization transition here reported are quite general for extended 
dynamical systems, since analogous behaviours have been recently observed for two coupled 
CMLs without any external noise EH . 



Moreover, as far as the control of chaos is concerned when the erratic behaviours present 
in the extended system are due solely to non-linear mechanisms (as it happens when the 
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maximal Lyapunov is negative, but V is still positive) the control schemes based on linear 
analysis [26| should fail and new "non- linear" methods have to be introduced. We believe 
that the appropriate indicator to employ in this context is the Finite Size Lyapunov Exponent 



27J , because it is able to capture infinitesimal as well as finite amplitude perturbation 
growth. 

We also expect that the observed phenomenology is not peculiar of CMLs and the present 
study can apply to a wider class of spatially extended dynamical systems, like coupled 
oscillator or PDE's. 
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FIGURES 



f(x) 




-0.5 



4 6 



FIG. 1. The map © for a = 4, b = 2 and c = 0.5. 
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FIG. 2. n versus L for coupled maps (||) for a = 4, 6 = 2, c = 0.5, with e = 1/3: a) exponential 
scaling observed for 2.49 < a < 2 A3; b) logarithmic scaling for a = 2.6 (filled circles). The values 
of r\ have been computed with A = 10 -8 — 10 -12 and averaged over 1000 - 5000 initial conditions. 
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FIG. 3. Behaviour of the propagation velocity V (circles) and of the maximum Lyapunov 
exponent A (triangles) versus noise amplitude a for coupled logistic maps with e = 2/3. Both 
quantities have been computed for L = 1024, averaging over 10 3 initial conditions each one followed 
for 0(1O 5 ) time steps. The dashed lines are a guide for the eyes. 
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FIG. 4. Exponential scaling of t\ versus L (reported in a lin-log scale) for coupled logistic 
maps with e = 2/3 at different values of the noise amplitude: a = 0.3 (filled circles), 0.32 (empty 
circles), 0.35 (filled diamonds), 0.38 (empty diamonds). The inset shows the logarithmic scaling of 
Ti versus L (in log-lin scale) for a = 0.45 . The values of t\ have been computed with A = 10~ 10 
and averaged over 10 3 initial conditions. 
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FIG. 5. (a) Exponential scaling of t\ versus L (reported in lin-log scale) for coupled maps (|5|) 
for a = 4, b = 1000, c = 0.5, with e = 1/3 and a = 2.1 (filled circles), (b) Logarithmic scaling of t\ 
versus L (shown in log-lin scale) for a = 2.5 (filled circles). In this case <ta = 1-99 and ay = 2.15. 
The values of t\ have been computed with A = 10~ 8 and averaged over 10 3 initial conditions. 




FIG. 6. (a) The inverse of the scaling factor l/£ as a function of a^ 1 — a^ (<ta = 2.5015) is 
reported for coupled (||) maps with b = 2.0 and e = 1/3 and 2.49 < a < 2.43; (b) l/£ as a function 
of a~ l — ay 1 {ay = 0.4018) for coupled logistic maps with e = 2/3 and 0.30 < a < 0.38. 
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FIG. 7. The spatio-temporal averages of the indicators £ (a) and p (b) are reported as a function 
of the noise amplitude a for for coupled logistic maps with e = 2/3. The data have been obtained 
by averaging over a time t = 10 and 20 different initial conditions, for L = 4096 and considering 
a threshold A = 10~ 12 . 




FIG. 8. The first passage time n is reported in a log-log scale as a function of system size L for 
coupled logistic maps with e = 2/3 and for various a. The data have been obtained by averaging 
over 3,000 - 25,000 different initial conditions and considering a threshold A = 10~ 12 . 
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FIG. 9. Logarithm in base 10 of the density of active sites as a function of log 10 (t) for coupled 
logistic maps with e = 2/3 and a = 0.4018. The data have been obtained by averaging over 10, 000 
different initial conditions and considering a threshold A = 10~ 15 . 
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FIG. 10. The number of "infected" sites N(t) is reported in a log-log scale as a function of time 
t for coupled maps d) with a = 4,6 = 2,c = 0.5 and e = 1/3 and for various a. The numerical 
data (filled circles) have been obtained for a chain length L = 20,000 by averaging over 1,000 
different initial conditions and considering a threshold A = 10~ 15 . The dashed line indicates a 
scaling of the Edwards- Wilkinson type (with exponent 1/z = 0.5) while the dot-dashed line has a 
slope 1/z = 2/3 analogous to that of the KPZ scaling. 





21 



